Exactly solvable model for isospin S=3/2 fermionic atoms on an optical lattice 
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We propose an exact solution of a model describing a low energy behavior of cold isospin S=3/2 
fermionic atoms on a one-dimensional optical lattice. Depending on the band filling the effective 
field theory has a form of a deformed Gross-Neveu model with either 0(7) x Z2 (half filling) or 
U{1) X 0(5) X Z2 symmetry. 



Although high symmetries do not occur frequently in 
nature, they deserve attention since every new symmetry 
brings with itself a possibility of new physics. According 
to 111, the SO (5) and SO (7) symmetries can be realized 
without fine-tuning of parameters in lattice models of 
isospin-3/2 fermions with a contact interaction. An ex- 
perimental realization was found in , where a four com- 
ponent gas of fermionic K^*^ atoms was kept in an optical 
trap for 50 seconds. In this letter we will limit ourselves 
to one dimensions where the corresponding models hap- 
pen to be exactly solvable. The appropriate model is the 
one-band generalized Hubbard model 

H = -tY^[i^+{^ + l)^b^{i) + h.c.]-^lY,i^ii')M^) 

i,a i,a 

+ U0 J2 P0,0i^P0,0{^) + f/2 ^ P^,mi^)P2,m{l) (1) 
i i,m 

where a and m take values a — ±3/2, ±1/2 and m = 
0, ±1, ±2. The singlet (total spin St = 0) and quintet ( 
St — 2) pairing operators Pqq, ^-^^ defined through 
the Clebsch-Gordon coefficients 

Q,/3 

where F = 0, 2 and m — —F, ~F + 1, ...F. In 1] it was 
shown that this model, together with its continuous ver- 
sion, have a precise U(l)xS0(5) symmetry at arbitrary 
filling and SO (7) symmetry at half filling. 

It was demonstrated later that, at least in the low en- 
ergy limit, the total symmetry of these models is greater 
being U(l)xO(5)xZ2 and SO(7)xZ2 respectively 110. 
Following the standard procedure (see, for example, p), 
at temperatures much smaller than the Fermi energy one 
can linearize the fermionic spectrum in the vicinity of the 
right and left Fermi points and to arrive to the effective 
field theory description. If the band filling is incommen- 
surate, the charge sector decouples from the other degrees 
of freedom at low energies. In this limit the Hamiltonian 
can be decomposed into a sum of the charge and the spin 
part H K, H charge + -ffspin, whcrc the charge sector is de- 
scribed by the Gaussian model. All nontrivial physics is 
concentrated in the spin Hamiltonian. It has SO(5)xZ2 



symmetry. The four-fermion interaction can be written 
as a sum of products of the SUi (4) currents of the right 
and left chirality though its symmetry is lower than 
SU(4). To make the true symmetry manifest one can 
employ the nonlinear transformation of fields suggested 
in 7] and rewrite the Hamiltonian density in the spin 
sector in terms of six species of Majorana fermions (the 
bosonized form of this Hamiltonian was obtained in 4]): 

i ^ 

'H = -l^^VsiXadxXa - XadxXa) + 

g (^^^XaXaj +25x(X0X0) (^^^XaXa^ ■ (2) 

where 5 ^ —{UQ + U2), gx ^ ~3U2 + Uo. One can run few 
simple checks to justify that above model does represent 
the continuum limit of QJ. First, the interaction in Q 
has the right scaling dimension (in other words, this is 
a model with four-fermion interaction). Second, at tem- 
peratures much higher than the spectral gap but much 
smaller than the bandwidth, when the interactions can 
be neglected, the two models have the same specific heat. 
This corresponds to model Q having the right ultraviolet 
central charge Cjjv = 3. Third, model 101 has the sym- 
metry predicted in (the Z2 symmetry corresponds to 
Xo —Xo)- In the limit g = gx the 0(5) XZ2 model be- 
comes the SO(6)=SU(4) Gross-Neveu (GN) model. This 
limit is also contained in original model J^l At half 
filling, due to the Umklapp processes, the charge mode 
is not decoupled. Therefore the entire effective Hamilto- 
nian is given by Eq.Q with seven degenerate Majorana 
fermions instead of five (the 0(7)xZ2 model). 

The 0(5) xZs (Eq. El and 0(7) XZ2 models (further 
down we refer to them as (5+1) and (7±1) models) have 
nontrivial dynamics for g > 0, when they scale to strong 
coupling and spectral gaps are generated. As field the- 
ories they exist in two limits, both corresponding to in- 
tegrable models 8]. One limit, associated to the largest 
symmetry (S0(6) or S0(8)), is reahzed when g = gx- 
The corresponding exactly solvable model is the S0(2N) 
{N ^ 3,4) GN model |lIll[ilE3- 

In this letter we discuss the exact solution of the 
anisotropic version of the models with g =/= gx- Most 



2 



of the calculations are conducted for the (5+1) 
the end of the paper we briefly discuss the 0(7) XZ2 gen- 
eralization. We construct the exact solution using the 
bootstrap procedure as described, for instance, in 
Namely, guided by the symmetry of the model and the 
perturbation theory results we suggest the two-particle S- 
matrix. In integrable models, where multi-particle colli- 
sions are representable as a sequence of independent two- 
particle ones (factorizability) , such S-matrix contains all 
information about the spectrum and the thermodynam- 
ics, as well as off-shell properties, such as correlation 
functions. The two-particle S-matrix must satisfy the 
crossing, unitarity and the Yang-Baxter conditions 
The latter one is a condition of associativity of the al- 
gebra of creation (annihilation) operators. This condi- 
tion, being an overcomplete system of equations for the 
S-matrix elements, is very restrictive. Ones such solu- 
tion is found, one has to complete the bootstrap process 
by calculating the free energy and comparing it with the 
perturbation theory for the model at hand. 

As we shall see, the bootstrap solution we obtain does 
not contain the fully symmetric (0(6) or 0(8)) limit. 
This suggests that the symmetric limit is unstable and 
cannot be reached asymptotically from the exact solution 
of the asymmetric model. This conclusion is supported 
by the analysis of the renormalization group (RG) equa- 
tions. The RG equations for the 0(2N-|-l)xZ2 model are 



gx = -2Ngxg. 



(3) 



They have the following RG invariant: C = [g^ ~ 
9x)/\9x\^,h = {2N - l)/N. From this it follows 
that the distance from the symmetric line \g — gx\ ~ 
C||5x|'^ "'^■'^^ does not decrease for all > 1 under 
the increase in gx rendering this line unstable. 

The integrability of the 0(2N+l)xZ2 model is a nat- 
ural extension of the = 1 case studied in detail in 
Refs. 0, In that case the two-particle S-matrix has 
the form 



3+1 — 



-1 



(4) 



where a, a are in and out spinon indices, (^{6) — (e^ — 
i)/(e^ -t_i) and £'^^(2) is the S-matrix of SU(2) Thirring 
model [ll[ [iTf . An important check of the validity of this 
S-matrix is that the Thermodynamic Bethe ansatz equa- 
tions constructed from it yield the correct UV central 
charge C3+1 = 2 and particle multiplicity. 

The S-matrix for the model was outlined in as 
a generalization of the (3+1) case. In this letter we give 
further details of this solution. The suggested S-matrix 



has the following form: 



S5+1 — 







-i 




Sss 


SsO 


-i 


SsO 


-1 



(5) 



The indices s and v label massive kinks belonging to the 
4-dimensional spinor representation and vector particles 
(associated to the original fermions Xa ( a = 1, ...5)). The 
kink's mass is M^; the vector particles have masses My. 
In the present model they are not bound states of kinks. 
There is an additional singlet particle (corresponding to 
the original xo fermion) with mass Mq. All masses are 
expenentially small in the bare coupling constants. Since 
there are only two couplings in the theory, Mg^My and 
Mq are not independent, though to determine this de- 
pendence one needs a microscopic derivation. As one 
might have expected, the Majorana fermion xo acquires 
a nontrivial phase factor by scattering on a kink: 



5,o(0) = ?(e)<5^ m 



■36 



+ i 



(6) 



(a, a are in and out kink indices) and 
matrix with the other vector particles, 
notations are as follows. S'"", iS"* and 
matrices of vector and spinor particles 
written as a sum of projectors. Pa, that 
product of two representations onto the 
resentation labeled a. For instance, the 
has the form 



has a trivial S- 
The rest of the 
S"' are 0(5) S- 
. They can be 
map the tensor 
irreducible rep- 
spinor S-matrix 



Pa. 



syrti 



6l + i7r/3 



Pv 



ITT 



(7) 



where Po,Pv,Pasym represent projectors onto singlet, 
vector and antisymmetric tensor representations. The 
form of the prefactors, 5*^^, Sgv, S^, can be inferred from 
the kernels (|15l) via 



Sabi^) = exp 



duj 



(8) 



and for the vector particles can be found in Ref.sl^ HH 
For the spinor particles we have 

Sss{e)^f{e)/f{-e), f{e)= (9) 

r(2 + a/c + 3i6l/27r)r(4 + 3fc + 3i6l/27r) 



n 

fc=0 



r(5/2 + 3fc + 3i6l/27r)r(7/2 + 3fc + 3i6i/27r) 



It has no poles on the physical sheet, thus, as we have 
already mentioned, spinors do not create bound states. 

The Bethe ansatz equations associated with the S- 
matrix ((SJ are 
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no 



&=1 



U, 



-iM^Lsinh^'j"'' _ 



b^l c^l p^l 

Us rflu 



Us 



6=1 



b=l 



-lAfni sin 



(lOb) 
(lOc) 



b=l 



where the numbers Us^v > ttIs^v and uq correspond to 



physical particles (spinors, 0(5) vector and singlet Ma- ing the standard methods we obtain the thermodynamic 
jorana fermions respectively) and vanish in the ground Bethe ansatz (TBA) equations for the (5+1) model. The 
state. The functions e„ have the form e„(a;) ~ {x + 
in7r/3)/(a; — in7r/3) and the auxiliary particles satisfy the 
eigenvalue equations for the 0(5) transfer matrix |l2| 



n ei(Aa - Mb) n ei/2(^a - Oi''>) = U ei(Aa - A,) 

6=1 c=l 6=1 

n ei(A^a - Afc) n ei(Ma - 4"^) - n ^^(Ma - Mc).(n) 



Taking the continuum limit of Eq.s (|10llll() and foUow- 
ig the standa 
ethe ansatz ( 
free energy is 

F/L^-T J2 — /d0cosh01n[l + e-'"W/^] (12) 



a—0.s,v 



where the dressed energies (a = 0, v, s) satisfy the 
nonlinear integral equations 



6=1 



6=1 



6=1 



and 



eo{0) = Mo cosh 6* - Ts^^/^) ^ j^j^ ^ g-e,(e)/Tj 

es{e) = M, cosh 61 + TY,s * ln[l + e-<^=(^)/^] + TY^y * ln[l + c-^-^^)/"^] - Ts^^/^) * ln[l + e-^°(^)/^] 
Tall/2) *ln[l + e-'"'/T] 

e^(0) = M„ cosh 61 + TYy^ * ln[l + e-'^^'^'>/^] + TY^^ * ln[l + e'^^f^)/^] - ra„ * ln[l + e""^- '/^] 



ln[l + e^i^^VT] _ ^(1/2) ^ _^ g-e(^)/T] _^ ^(1/2) ^ , J^f^ ^_ g-.i^'/T] ^ ^(1/2) ^ j^f^ ^ ^-eJ7 

ln[l + c^" '/^] - A„,,„ * ln[l + e-^"'/^] + * A„,™/2 * ln[l + e"'™'/^] = a„ * ln[l + c^^"/^] 

I 



(13a) 

(13b) 
(13c) 

(14a) 
(14b) 



These equations can be solved numerically to get a de- J dii exp 
tailed description of the thermodynamics. For the pur- 
poses of this paper it will be sufhcient to study asymp- 
totics of the free energy at large and small tempera- 
tures. In the above equations * stands for the oper- 
ation of convolution / * g{9) = J duf{9 — u)g(u) and 
all kernels are related to the Fourier images as f{d) = 



Ysv (w) 



. The kernels in Fourier space are 

e-2|-l 



4cosh(3w/2)cosh(u;/2) 

el'^l 



(15) 



~ cosh(u;/2)el 
Yyy{uj) = tt:: — — 1 



cosh(3w/2) 

cosh((jj/2) 



coshc 



2cosh(3w/2) 
a„{uj) = e""l'^l, s{lu) = l/2coshw, ^/^(w) = 

^rim(w) = COth |w| 

and we used the notation 

/(i/2)(^)^/(^/2). (16) 

At Mo = {gx — 0) model Q decouples into one mass- 
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less Majorana mode and the 0(5) GN model. As usual, 
we expect that the massless mode does not appear in 
the Bethe ansatz for massive particles. We observe that 
in this limit the above TBA equations coincide with the 
ones derived in ^ for the 0(2N+1) GN model for the 
case N = 2. The dressed energies cq, Cy, £« correspond 
to eN,-i, eAf~i,0) ^N,o of This is an important self- 
consistency check of our solution. Please note that while 
eat.-i was an auxiliary particle in the 0(5) GN model, in 
our model {Mq ^ 0) it is a physical particle. 

With the TBA equations one can run two additional 
checks. First, by calculating the free e nerg y at T >> Ma 
(for this end one can use the results of |l^) one finds the 
universal asymptotics F/L = — Ct/yT^7r^/6 such that at 
Mq ^ one gets Cj/y = 3. This is the correct asymp- 
totics for a theory of 6 Majorana fermions. On the other 
hand the T « Ma limit reproduces another asymptotics 

F/L = -TY^TaMa I ^cosh0e-*^"^°^h9/T 

with = 4, r„ = 5, ro = 1 being the correct particle 
multiplicities. 

The generalization for the (7-1-1) model is straightfor- 
ward. Here we have 8 kinks, 7-1-1 vector particles and 
28 tensor particles. The S-matrix has the same block 
structure as (jSJ, but with 

^'^^^ ^ ^ 

and an additional block Su for the tensor particles. The 
blocks Syy and Su coincide with the vector and tensor S- 
matrices of the 0(7) GN model and are described in [l^ . 
The spinor block Sss is similar in its form to Q with 
7r/3 substituted by 7r/5. It does not have poles since the 
vector particles are not bound states of spinors. 

The spin sector at arbitrary filling and the entire model 
at half filling are quantum liquids. Away from half 
filling the corresponding order parameters and the phase 
diagram are described in P[ At moderate forward 
scattering there are two phases corresponding to Density 
Wave or the superfluid state of the BCS type. They 
are distingushed by the sign of coupling gx which does 
not affect the excitation spectrum. Both these phases 
have power law correlations and are separated by the 
Ising type Quantum Phase transition. Generalizing the 
arguments given in 0|, we conclude that for half filling 
there are also two phases. Just one of them has an order 
parameter local in terms of the fermions: it is 2fcF = tt 
Charge Density Wave (CDW). It condenses at gx < 0; 
at ffx > the order parameter is nonlocal corresponding 
to a hidden (topological) order. 

The fact that the 0(6)( or 0(8)) symmetry is broken 
down to 0(5)xZ2 (or 0(7)xZ2) will play a role for the 
correlation functions. Indeed, the 0(6) and 0(8) groups 
have two irreducible spinor representations and 0(5) and 



0(7) have just one. As it was shown in |l9|,in theO(2N) 
GN model right and left moving fermions transform ac- 
cording to different spinor representations making it im- 
possible to have a nonzero Green's function containing 
right and left fermions. However, since the 0(8) symmet- 
ric point in model ^ is unstable, the number of kinks 
in the exact solution is the same as for the 0(7) group, 
which has just one irreducible representation. Therefore 
the Green's function of the right and left fermions may 
be non zero. It is logical to assume that it does exist in 
the dimerized phase corresponding to gx < 0. We intend 
to discuss this issue in subsequent publications. 
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